Critical properties of doped coupled spin-Peierls chains 
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Using numerical Real Space Renormalisation Group methods as well as Stochastic Series Expan- 
sions Quantum Monte Carlo simulations a generic model of diluted spin-| impurities interacting 
at long distances is investigated. Such a model gives a generic description of coupled dimerized 
spin-Peierls chains doped with non-magnetic impurities at temperatures lower than the spin gap. 
A scaling regime with temperature power-law behaviors in several quantities like the uniform or 
staggered susceptibilities is identified and interpreted in terms of large clusters of correlated spins. 
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Low dimensional gapped quantum magnets have at- 
tracted a lot of interest in condensed matter physics for 
many years. The possibility of doping such systems has 
lead to an extremely rich emerging field. The discovery 
of the first non-organic spin-Peierls compound CuGcC>3 

and, soon later, its doping with static non-magnetic 
impurities realized by direct substitution of a small frac- 
tion of copper atoms by zinc or magnesium Q atoms 
offered a new challenge for the theorist and an ideal ex- 
perimental system to test new theoretical concepts. 

Our aim here is to analyse the low temperature prop- 
erties of a typical two-dimensional array of coupled, 
frustrated and dimerized antifcrromagnetic (AF) spin- 

1 chains doped with non-magnetic impurities. For that 
purpose we use the low-tempcraturc effective model de- 
rived in previous work: each non-magnetic dopant re- 
leases a spin-i, localized in its vicinity These effec- 
tive spins become the only remaining low-energy degrees 
of freedom at temperatures lower than the spin gap i.e. 
the energy scale of condensation of the background spins 
into nearest-neighbor dimers. The corresponding effec- 
tive Hamiltonian describes interacting spins \ randomly 
distributed on a square lattice (of size L x L), 
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where the occupation number e r takes random values 1 
(0) with probability x (1 — x), x = N s /L 2 is the dopant 
concentration, N s the number of spins \. The effective 
interaction J(r) computed by Lanczos exact diagonaliza- 
tions from the original microscopic model Q bears im- 
portant properties: (i) it has opposite signs on the two 
(relative) sublattices and, hence, is non-frustrating in na- 
ture and (ii) it displays a typical (spatially anisotropic) 
exponential (long distance) behavior, characterized es- 
sentially by two length scales, £m and 
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where Ax = Ax + x max (see later) and Ax = \x\ — x%\ 
{Ay = | j/i — 2/2 1) is the separation between 2 dopants in 



the longitudinal (transverse) direction. Typically, we as- 
sume hereafter £|| = 2.5 and £i_ = 1. Note that Eq.JJJ) is 
generic for most doped spin gapped systems 0, In 
order to use a consistent and realistic description of the 
actual experimental compound, we shall use also the spe- 
cific form derived in Ref. Q at short distances for an AF 
interaction i.e. when Aa; + Ay is odd J {Ax, Ay) oc 
Ax up to a maximum at distance x max = 2£|| (x max = 
for Aa; + Ay even). At larger distances, the asymptotic 
form of Eq.J2J) is used. 

As supported by experiments 0, the present frame- 
work implies naturally that the low-temperature uniform 
susceptibility scales with the impurity concentration and 
displays a Curie- like behavior x = C/T. As argued in 0, 
due to the presence of AF as well as ferromagnetic (F) 
couplings [lOj, large (weakly coupled) clusters of spins 
form. Using a simple classical random walk argument, 
the effective spin S eS of the clusters containing n spins 
is simply given, on average, by {S cS ) = ^(n) 1 / 2 . Con- 
sequently, the Curie constant (per spin) should saturate 
to a value equal to i((5 cff )) 2 /(n) = 1/12 when T -> 0. 
Preliminary numerical results are in agreement with this 
prediction Here we shall use both numerical Real 

Space Renormalization Group (RSRG) [T3. Il3j supple- 
mented by Stochastic Series Expansion (SSE) Quantum 
Monte Carlo (QMC) Q to investigate the critical prop- 
erties as one approaches the zero temperature limit. 

For a single spin-^ chain with random AF couplings, 
Fisher demonstrated the existence of a universal fixed 
point for the renormalization group transformation |l 3t | , 
thus providing a strong argument in favor of such a pro- 
cedure. In the case of both random F and AF couplings, 
while there is no exact solution, the existence of a fixed 
point was shown unambiguously by Westerberg et al. flflj | 
using a RSRG approach. Its validity was discussed very 
carefully in 01 an d the predictions were successfully 
checked by QMC 0|. In the similar 2D case we ad- 
dress here, the distribution of couplings P{\ J\) behaves 
like 1/| J | and does not depend on the dilution x (apart 
from logarithmic factors and a finite size cutoff J m i n ). 
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The RSRG is expected to work particularly well for such 
a singular distribution. 
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FIG. 1: Schematic picture of the doped spin-Peierls system. 
Thick bonds stand for dimers and non-magnetic impurities 
(released spins |) are represented by open circles (black ar- 
rows). The initial RSRG step is illustrated starting from a 
typical configuration with 4 impurities: (a) The strongest cou- 
pled pair is identified (red line), (b) This pair, e.g. ferromag- 
netic here, is replaced by a spin S = 1 (red arrow), (c) The 
couplings with all other spins (dashed lines) are renormalized. 



Following the pioneering work of Bhatt and Lee (l2^ . 
we extend the RSRG scheme to hamiltonian with 
F and AF long-distance couplings 0, Let us de- 

fine the effective interaction as Jj j where i and j la- 
bel the randomly distributed spins and run from 1 to 
N s . One single RG step is described as follows : 1) 
Identify the most strongly coupled pair of spins (Si, Si) 
i.e. with the largest energy gap Ai i2 , A ii2 = ^1,2(1 + 
15*1-5*21) if Ji, 2 > (AF) and A h2 = -Ji, 2 (Si + 
S2) if Ji,2 < (F). Note that Ai j2 defines the energy 
scale of the transformation. 2) Replace it by an effective 
spin 5 = I Si — S2 1 if the coupling is AF or 5 = Si + S2 
in the F case. 3) Renormalize all the magnetic couplings 
with the following rules : (i) If 5 7^ 0, as given by a first 
order perturbation theory, the new couplings between S 
and all the other spins (53, 54, Sn 3 ) are set to 
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(ii) If 5 =0, the pai r (5i,5 2 ) is frozen. Using a clus- 
ter approximation [12| that involves only the extra pair 
(53, 54) the most strongly coupled to 5i and 5 2 and a 
second order perturbation, the coupling J3 4 is renormal- 
ized as 
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(Ji,3 — ^2,3X^2,4 — Ji,i)- 



The same procedure is then reiterated. We also check 
that the RSRG preserves the non-frustrated character of 
the problem. 

Due to the presence of both F and AF couplings, clus- 
ters with large effective spins are created during the 
procedure similarly to what occurs in the ID random 



F-AF spin-i chain [ljj. At each RG step, the en- 
ergy scale A decreases and both the number of inac- 
tive spins frozen into singlets and the number of clus- 
ters build from a large number n correlated spins- \, in- 
crease. The aforementioned random walk argument pre- 
dicts that, the average number (n) of spins- ^ inside clus- 
ters and their average spin mag nitudc (5 cff ) should be 
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)*/ 2 at low enough temperatures. 



Therefore we expect the effective spin of these clusters to 
grow monotonously as the energy scales down. We have 
analyzed this process using the RSRG scheme to com- 
pute both (5 ) and (n) as a function of (Aq). This is 
shown in Fig.|21which clearly demonstrates the formation 
of large moments. Moreover, power-law divergences like 
(5 cfr ) - (Aq)-^ and (n) ~ (A )~ k{x) are observed 
with k — lot. We have plotted the behavior of these 
exponents vs x in Fig. [7] Interestingly enough, we find 
that a depends on x in contrast to the random F-AF spin 
chain for which a = 0.22±0.01 llfl. Note however that 
the universality class identified in [lfj as a RG fixed point 
occurs only for initial gap distributions less singular than 
P C (A) ~ A~ Vc , with 0.65 < y c < 0.75. For more singular 
distributions, as it is the case here, critical exponents are 
not universal anymore, at least in ID [i"fl |. 




FIG. 2: Average effective spin (S eS } of the clusters of active 
spins vs the energy scale (Ao) for six different concentrations 
x indicated on the plot. Numerical RSRG data obtained for 
N s = 1024 spins over more than 10 4 samples. Inset: for 
the same samples, average number (n) of initial spins- \ per 
cluster vs (Ao). Dashed lines are power-law fits (see text). 



We now turn to a brief discussion about the validity of 
the RSRG approach w.r.t. the large (formally infinite) 
connectivity of the model. Whereas the moments eventu- 
ally order AF at zero temperature for arbitrary small con- 
centration 01 ? at finite temperatures the physics is dom- 
inated by short range couplings, i.e. inside a domain de- 
fined by Ax < £|| and Ay < Consequently, one could 
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define a "short-range connectivity number" z sr which is, 
for a given spin, the average number of neighboring spins 
belonging to such a short range couplings region. It is 
easy to check that z sr = Af sr x(l — N" 1 ) ~ AT sr x where 
M S r is the number of sites contained in the short-range 
region. For the present model AT sr — 65 sites as seen 
numerically. The evolution under the RSRG scheme de- 
pends on the value of z sr , hence on x. When, let us say, 
z sr < 1 (i.e. x < 1.55%), the probability to have extra 
spins strongly coupled to the most strongly coupled pair 
(Si,^) remains very small, so that the necessary con- 
dition J1.2 ^> J1.i1J2.il Vi > 2 is fulfilled. On the other 
hand, if z sr > 2 (i.e. x > 3.1%) a "percolation thresh- 
old" is reached and one expects the formation of strongly 
coupled clusters reaching the system size. 
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FIG. 3: Curie constant per spin (C) plotted vs the energy 
scale, (a) Quantum Monte Carlo SSE results shown vs T for 
N 3 — 256 spins and 9 different concentrations x indicated 
on the plot. Disorder is performed on 10' ! to 10 4 samples, 
(b) Numerical RSRG results shown for N s = 1024 spins and 
6 different concentrations x as indicated on the plot, vs the 
RG energy scale (Aq). Error bars are typically smaller than 
symbol sizes, the number of samples always exceeding 10 4 . 
The full line correspond to the saturation value of 1/12. 



Since the above arguments are still qualitative |l7j . 
we have confronted the RSRG results to QMC data ob- 
tained on the same model and for the same parameters. 
The SSE method 0] supplemented by the /3-doubling 
scheme 0] , already used in 01 > can reach extremely low 
temperatures. Figs.OHshow results for the Curie constant 
per spin obtained with both methods. The RSRG com- 
putation of C is performed, at each RG step, using the 
formula C = j^- N cr a(cr + 1) where N a is the num- 
ber of active effective spins of size a, the data being then 
averaged over disorder. We have also checked that finite 
size effects are negligible when N s > 256 and we have 
chosen N s = 1024 in most computations. We observe 
a qualitative agreement between Fig^a) and Fig|3jb) 
where, at high temperatures, the spins behave as param- 
agnetic free magnetic moments (giving a Curie constant 
of j per spin) and where saturation to ^ is observed at 
low T, the spins being correlated inside large clusters. At 



small concentration, the agreement becomes even quan- 
titative, as seen in the Fig0]where a comparison between 
RSRG and SSE is shown for the six lowests values of x. 
The agreement remains quite good for larger concentra- 



te 




10" 

T, <A ) 



FIG. 4: Direct comparisons between RSRG (full lines) and 
SSE simulations (symbols) of (C) — 1/12 are shown for the 
6 different concentrations indicated by (a), (b),..., (f) vs the 
RG energy scale (Aq) or the SSE temperature T. 



tion up to x ~ 1.56% corresponding to z sr ~ 1, hence 
corroborating the qualitative argument stated above. 

We now turn to the analysis of the scaling regime. At 
very low temperatures, the quantum corrections (C(T)} — 
1/12 are expected to behave like T 7 where 7 = a in 
ID 0. In the present case, the procedure used in Figs.|U 

T 7 sse 
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FIG. 5: Curie constant (C) plotted vs T 7 sse for the QMC 
data in (a) and vs (Ao) 7r g for the RSRG data in (b), in 
order to obtain the best data collapse at low energy. 



to extract the doping dependence of 7 from the QMC and 
RSRG susceptibility data is the following: First, one esti- 
mates the value of 7 at the lowest concentration available, 
i.e. iK mi n ^ 0.41% for the SSE and x min = 0.1% for the 
RG, via a direct reliable power-law fit, giving 7 SSE = 0.12 
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and 7 H 



0.065 in thoses cases. Then, the other esti- 



mates for larger x are determined in order to obtain the 
best collapse of all the data plotted vs T 7 on a universal 
low temperature curve, as shown in Figs.El We note that, 
for large x, deviations at high temperatures can be at- 
tributed to a transient and plateau regimes also visible in 
Fig. |3fa) for x > 1/36 |2(|. Similarly, one expects for the 




FIG. 6: T(xsta g (r)) plotted vs T for five different concentra- 
tions. Full lines are fits corresponding to power-law behaviors 
~ T~ 27 '. All data are computed by QMC using N s = 256 
random spins and averaged over disorder. 



staggered susceptibility (per spin), (XstagCUl) T -1 ~ 27 . 
In ID, 7' is expected to be equal to a |l9j, but in the 
present case, by direct fits of the low T (see Fig. [JjJ, 
we found 7' ~ 7 ~ 2a. Thoses exponents are plotted 
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FIG. 7: Exponents ct(x), 7(1) and 7' (a;) extracted from var- 
ious SSE and RSRG data as indicated on the plot (see text 
for details). Straight lines are ~ y/x. 



in Fig. [7| An overall very good agreement is seen be- 
tween the different methods, in particular between the 
estimates of 7 obtained from the analysis of the Curie 
constant computed by QMC and RSRG. We stress again 
that the exponent a deduced from the analysis of the 
change of cluster sizes and spins is roughly a factor of 2 
smaller than 7. Interestingly enough, 7 ~ 2a oc yfx in 
all cases. 

To conclude, confronting results obtained by both 
state-of-the-art QMC simulations and numerical RSRG 



methods, we have achieved a physical understanding of 
the critical properties of a generic 2D model of diluted 
S = I spins with long-ranged interactions. For small 
dopant concentration, the RSRG appears to be in ex- 
cellent agreement with QMC computations. Power-law 
temperature behaviors in several quantities like the uni- 
form or staggered susceptibilities are revealed and inter- 
preted in terms of large clusters of correlated spins which 
also have interesting scaling properties with temperature. 
It is remarkable that such properties are observed above 
the T = ordered magnetic groundstate. Moreover, it 
is the first example of a two-dimcnsionnal random mag- 
net exhibiting a large spin phase with a disorder (the 
concentration x) dependence of the critical exponents. 
Whereas doped CuGeOa might be a good candidate for 
the experimental observation of a critical regime, the so 
large three-dimensionnal ordering temperature prevents 
such an observation. More strongly diluted samples are 
necessary to reach the scaling regime. 
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